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C  -  absolute  concentration  difference  at  a  field  point 
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INTRODUCTION 

Recent  experiments  on  cavity  pressure  measurement  [2,3]  have 
demonstrated  that  the  cavity  pressure  inside  a  cavity  varies  with  cavity 
oscillations  and  that  the  cavity  pressure  fluctuations  are  significant  at 
lew  frequencies.-  The  gasecus  diffusion  across  a  cavity  surface  naturally 
influences  tho  cavity  pressure  and  hence  the  characteristics  of  the 
cavity  (viz.,  size,  shape  and  behavior)  to  a  great  extent.  While 
previous  investigations  have  addressed  the  steady  gaseous  diffusion  [1,4] 
across  cavity  surfaces,  a  corresponding  analysis  for  unsteady  diffusion 
in  an  unsteady  cavity  flow  is,  to  the  best  of  our  knowledge, 
non-existent ., 

Unsteady  gaseous  diffusion  into  a  cavity  may  result  from  two 
distinct  mechanisms.-  First,  since  the  cavity  length  fluctuates,  the 
gaseous  diffusion  into  the  cavity  fluctuates  correspondingly.  Second, 
since  the  cavity  gas  pressure  itself  fluctuates,  the  concentration 
gradient  across  tna  cavity  surface  is  a  function  of  time..  In  this 
manuscript,  we  formulate  and  solve  analytically,  the,  unsteady  gaseous 
diffusion  across  a  two-dimensional  unsteady  cavity  surface.  The  results 
of  this  analysis  will  give  a  batter  understanding  of  the  parameters  tKv- 
affect  the  unsteady  gaseous  diffusion  across  a  cavity  surface. 


Assumptions 
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It  is  assumed  in  the  foregoing  analysis  that  the  flow  is 
incompressible,  two-dimensional,  isotropic  and  turbulent..  [Turbu.ence  is 
inherently  three-dimensional  in  nature.  Nonetheless,  a  two-dimensional 
flow  field  with  unit  depth  may  be  considered  for  diffusion  analysis.] 
Soussinesq's  hypothesis  [4|  is  used  in  turbulence  mcdelin-  and 
G.  I.  Taylor's  [5]  statistical  theory  of  turbulence  is  used  in  estimating 
the  mass  diffusivity  for  the  flow  field.  Ic  is  assumed  that  the 
mechanism  for  turbulent  momentum  transfer  and  turbulent  concentration 
transfer  are  identical..  It  is  also  assumed  that  the  gaseous  diffusion 
occurs  instantaneously  across  the  cavity  surface  and  that  the 
concentration  gradient  on  the  cavity  surface  fluctuates  in  phase  with  the 
pressure  fluctuations  inside  the  cavity. 

Theoretical  Analysis 

Gaseous  diffusion  will  occur  when  there  exists  a  dissolved-gas 
concentration  difference  between  the  free  stream  ana  the  liquid  on  the 
cavity  surface.:  If  denotes  the  dissolved  gas  concent  in  the  free 
stream  in  parts  per  million,  by  moles,  then  by  Henry's  law,  the  maximum 
partial  pressure  cf  gas  in  the  free  stream  PpS  is 


PFS  -  V  • 


(U 


D 
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where  $  is  the  Henry's  law  constant..  If  PG(t)  is  the  instantaneous 
partial  pressure  of  non- condensable  gas  in  the  cavity,  then  the  mean 
concentration  difference  c(t)  expressed  in  moles  is 

PG 

c(t)  -  (aj  -  --)  ..  (2) 

We  assume,  as  is  customary,  that  the  mechanism  for  turbulent 
diffusion  of  gas  in  the  liquid  and  across  the  wall  and  the  turbulent 
momentum  transfer  are  similar  and  that  gradients  in  mean  velocity  can 
be  neglected.  Then,  the  turbulent  diffusion  is  due  solely  to  the 
gradient  of  the  mean  concentration  c(x.y}..  Consider  a  turbulent  flow 
field  having  a  uniform  mean  velocity  Ua,  in  the  positive  X  direction, 
as  shewn  in  Figure  l.:  Let  the  instantaneous  concentration  at  a  point 
(x,y)  be  C(x,y,t)..  Also,  let  Cm  be  the  concentration  at  large  distances 
from  the  origin  and  0o  be  the  saturation  concentration  on  the  cavity 
wall.;  Let  c  -  C  -  Cgj  be  the  concentration  difference..  From  Fick's  law 
of  diffusion  and  conservation  of  mass  (see,  for  example,  Batchelor  [6]) 
it  follows  that 

(U  «  V)c  +  &  -  «V2c  (3) 

where  denotes  the  Cartesian  Laplaciar.  operator  and  «.  is  the  molecular 


diffusivity . 
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Let  the  instantaneous  velocity  components  be  represented  by 
(!!„  +  u'),  v'  and  w'  in  the  three  orthogonal  directions  and  let 
c  -  c  +  c'  where  c  is  the  time  mean  concentration  difference..  The 
components  u' ,  v'  and  w'  are  the  instantaneous  velocity  fluctuations 
and  c'  is  the  instantaneous  concentration  difference.  Then,  Eq.;  (3) 
becomes  (see,  for  example,  Goldstein  [7]) 


ic 

at 


+  u 


i£ 

®  dx 


2- 
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Time -averaging  the  above  equation,  we  obtain 


i£ 

at 


T  & 

<*>  dx 


i_ 

ax 


u'c') 


i_ 

ay 


u'c') 


w’c') 

(5) 


where  the  overbar  indicates  time  averaging.  Equation  (5)  shows  that  in 
turbulent  flows  the  eddy  concentration  transport  terms  -  u'c',  -  v'c'  and 
-  w'c'  add  to  the  molecular  concentration  transports  k  k  and 

k  respectively..  The  eddy  transport  terms  are  large  compared  to 
molecular  transport  terms  (7).  Ther.  fore,  the  latter  are  neglected  and 


Equation  (5)  reduces  to 
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i i 

eft 


+  u 


is  _  .  £_ 

®  3x  3x 


(U'c') 


oy 


(-  u'c') 


3z 


(-  w'c') 


(6) 


We  now  restrict  our  study  to  tvo- dimensional  flows..  Ve  represent 
the  turbulent  shear  stress  pu'c'  and  p u'c'  using  Boussinesq's  hypothesis 
as 


-  pu'c' 


(7a) 


and 


\  § 


(7b) 


where  A-j-  is  the  exchange  coefficient  or  eddy  mass  conductivity.  Let  u 
denote  the  eddy  mass  diffusivity.  Then, 


»  -  at/p 


(8) 


where  p  is  the  density  of  the  liquid.;  Substitution  of  Equations  (7)  and 
(S'1  in  (6)  yields 


1  oc  _o  £c  fils.  d£_ 

v  5  c  w  3x  “  ,  2  .  ? 

3x  oy 


(9) 


We  normalize  the  distances  along  x  and  y  directions  by  a 
characteristic  length  b  and  the  time  by  a  reference  time  b/U®- 


such  that 


t 

(b/U  ) 
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(10a) 

(10b) 

(10c) 


Then,  Equation  (9)  may  be  written  as 


.  £c  1_  slk  sic 

3t  ax  2R  .  2  a  2 

ax  ay 


(U) 


where  the  quantity  R  is  a  nondimens ional  number  given  by 


R  — 


U  b 

CO 

2u 


(12) 


and  the  quantities  x,  y  and  t  are  now  nondimens  ional 

The  boundary  conditions  pertinent  to  the  two-dimensional  unsteady 
diffusion  are: 


|p-0  y-0  ,  x  <  0  ,x>  Z(t) 

and 


(13a) 


c  -  c  (t)  ;  y-0  ,  0  <  x  <  £(t) 

o 


(13b) 


where  JJ,(t)  is  the  unsteady  cavity  length  normalized  with  respect  to  the 
characteristic  length  b  and  is  given  by 

H(t)  -  2.  +  t  sin(kt  +  <(>) 

o 


(14) 
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In  Equation  (14),  20  is  the  steady  cavity  length.;  c  is  the  amplitude 
of  :avitation  ^ illation,  «  -  «(k ,a)  and  k  is  the  reduced  frequency 
defined  by 

k  -  (2irf)bA'  .  (15) 

co 


v'r.e.c  f  is  the  frequency  of  body  oscillation  in  Hertz.  $  is  the  phase 
angle  between  the  body  motion  and  the  cavity  motion,  <p  “  $(k)„  The 
solution  of  Equation  (11)  subject  to  the  prescribed  boundary  conditions 
1,13)  is  obtained  in  stages.  First,  a  point  source  solutioi  that  has 
the  proper  radial  symmetry  and  behavior  is  obtained.  Second,  the  p^int 
source  solution  is  used  to  formulate  an  integral  equation  relating  the 
concentration  at  a  field  point  and  the  mass  flux  per  unit  length  along 
the  cavity.  Third,  the  concentration  difference  at  the  cavity  surface 
is  used  to  solve  the  integral  equation. 

The  fundamental  solution  to  Equation  (11)  for  an  oscillating  unit 
source  at  the  origin  is  [see  Appendix  C-l] 


c  ( x ,  y ,  t ) 


1 

2jtAt 


2,  Rx  i(kt+<£-wx) 
y  )  e  e J 


(16) 


Let  f(x,C)  be  the  instantaneous  mass  flux  per  unit  length  between  x  -  0 
and  x  -  11(c)  The  the  function  c(x,y,t)  can  be  expressed  in  terms  of 
f ( x ,; t )  as  follows 


c(x,y. t) 


eii. 

2s-A^ 


*(t) 


f(£.t)eR(x‘rKo(R 


(x  -  O 


2 

+  v  )e 


jk(t-x+Odf 


0 


(17) 
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If  we  now  require  that  c(x,0,t)  -  c  (t)  in  the  interval  0  <  x  <  x.(t), 

o 

the  source  strength  is  determined  from  the  integral  equation 


c 

o 


fij (kt+^) 
2^ 


[/  f($,t)eR(x*€)Ko(R(x  -  C))eJk(X‘0dC 
0 


i(t) 

+  f  f(S,t)e‘R(f‘X)Ko(R(£  -  x))ejk(^‘x)d?]  .  (18) 

x 

where  the  positive  branch  of  the  square  root  has  been  taken..  The 
function  c(x,y,t)  also  satisfies  the  conditions  that  c  -  0  as 
(x^  +  y^)  ®  and  that  “  0  when  x  <  0  and  x  >  £.(t). 

In  order  to  determine  f(x,t),  we  shall  make  use  of  the  fact  that 
R  is  a  large  number  for  the  present  study.;  Therefore,  K0  in 
Equation  (18)  can  be  replaced  by  the  first  term  of  its  asymptotic 
expansion, 

Ko(z?  -  e'2VfF  -  09) 


A  further  simplification  can  be  obtained  in  the  second  integral  in 
Equation  (18)  by  noting  that  the  strong  negative  exponential  will 
cause  f(£,t)  to  be  a  slowly  varying  function  compared  to  the  term 
e  R<“  '  ‘v;Ko(k(£  -  x)),  so  that  f(£,  t)  contributes  to  the  integration 
inly  near  £  -  x..  Therefore,  we  replace  tie  second  integral  in 
Equation  (18)  by  the  approximate  value, 


*jmj  y  u* »  J/* j  ■gqypwwwwwwwwww  w  i  jj  m  n  n  w*iwiww*t 
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?r  Jk 

^  f(x.c)  _ : _ 


J 2R  -  jk 


With  these  approximations,  Equation  (18)  may  be  written  as 


2_  Pj(kt^~J  ,  .  f 

2R  2xA^. 


I*  _  «$ +  - j 

n  A  -  £  72R  -  jk 


(20) 


If  the  Laplace  Transform  of  f(x,t)  is  denoted  by  F(s,t),  we  can  transform 
Equation  [20)  to  find 


F(s,t)  - 


2V2R  coe-j(k£^> 


sf 


_L 


(s  +  jk)1/2  (2R  -  jk)1/2 


On  the  cavity  surface,  from  Equation  (C.1.2), 


c  (t)  -  c. (x,0)e^ (kc+<^ 
o  I 


Equation  (21)  may  be  rearranged  to  read 


F(s)  -  2A^  J 2R  cx  J 2R  -  ik  J. 


s+Jk 


(21) 


(22) 


s[72R  -  jk  +  Js  +  jk] 

The  inverse  of  this  transform  gives  [see  Appendic  C.2]  for  f(x)  the 
result 


f(x)  -  Q  [“  -  ~  e  aXerfcb7x  -  erfAxj 
a  a  a 


(23) 
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Q  -  2At  72R  cx  J2R  -  jk  (24b) 

a  -  jk  (24b) 

b  -  /2R  -  jk  (24c) 

a  -  (2 jk  -  2R)  (24d) 


and  erf(x)  and  erfc(x)  denote  the  error  function  and  complimentary  error 
function  respectively.  We  note  here  the  fact  that  the  function  f(x) 
which  represents  the  mass  flux  per  unit  length  is  independent  of  time. 

The  rate  at  which  the  mass  is  diffused  per  unit  widt .  along  the 
entire  length  of  the  permeable  interval  is  obtained  from 
Jl(t) 

~  -  /  f (x , t)dx  (23) 

0 

i  •  e  < , 

*(t)  2 

-  Q  /  (-  .  —  e  aXerfcb7x  -  erf./ax)dx  (26) 

at  c  a  a 

0 


The  integral  in  Equation  (26)  has  been  evaluated  in  Appendix  C.3  ar.d 
the  result  is 


erfcbyi  +  ^—z  erfja^  -  —  -  bjal  erfjal 


erfVaJ,  ] 


(27) 
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Equation  (27)  gives  the  instantaneous  mass  flow  rare  across  a  two- 
dimensional  ca,rity.  In  order  to  utilize  this  equation,  the  variation  of 
cavity  length  as  a  function  of  frequency  of  oscillation  must  be 
established  [see,  for  example,  Chapter  3  of  Ref..  (3)]..  In  the  limit, 
when  the  reduced  frequency  k  is  zero,  it  can  be  shown  that  the  expression 
for  mass  flow  rate  reduces  to  the  steady  state  mass  flow  rate  established 
by  Parkin  (1). 

Conclusions 

We  have  developed  an  analytical  model  for  the  prediction  of  gaseous 
diffusion  across  a  two-dimensional  unsteady  cavity  surface.  This  model 
takes  into  account  the  change  in  cavity  length  as  well  as  changes  in 
cavity  pressure  in  predictin'-,  the  gas  diffusion  across  the  cavity 
surface..  The  expressions  for  mass  flux  per  unit  length  and  mass 
diffusion  rate  across  the  entire  cavity  length  reduce  to  those  obtained 
by  Parkin  (8)  when  the  reduced  frequency  k  is  zero.  We  reiterate  here 
the  fact  that  this  unsteady  gaseous  diffusion  analysis  is  valid  only  for 
harmonic  variations  in  cavity  length  and  cavity  gas  pressure.. 
Nonetheless,,  this  analysis  can  easily  be  extended  to  encompass  a  general 
nonharmonic  motion  of  the  cavity  by  Fourier  representation  of  the  cavity 
length  and  cavity  gas  pressure  fluctuations. 


15 


31  July  198i 
KR:lhz 


APPENDIX  C.l 


Fundamental  solution  to  Equation  ''ll)  for  an  oscillating  unit  source 
at  the  origin: 

Equation  (4.11)  is 


£c  L  ,  slk  ,  $-.9 ' 

at  dx  2R  ,  2  .  2 

3x  3y 


(C.1.1) 


Since  (t)  is  periodic,  we  assume  c(x,y,t)  to  be  periodic  and  write 


c(x.y.t)  -  c1(x,y)e 


j  (kt  +<P) 


(C.l. 2) 


where  is  the  ohase  angle  between  the  motion  of  the  body  and  the  motion 
of  the  cavity  and  j  -  l.:  It  is  to  be  noted  that  only  the  real  part 

of  the  analysis  that  fellows  is  pertinent.  Substituting  Equation  (C.1,2) 
in  (C.1.1),  we  obtain 


3c 


2c.jkR  +  2R  r-*  -  V  c, 
1J  3x  1 


(C.l. 3) 


To  solve  Equation  (C.l. 3)  we  replace  the  Laplacian  on  the  right-hand 
side  by  its  counterpart  in  plane  polar  coordinates  when  variations  only 
in  the  radial  direction  are  permitted.  Then,, 


.2 

3c.  Q  c.  .  , 

2c.jkR  +  2Rri-  — r1  +  1  “ 
1J  3x  ,2  r  3r 

at 


(C.1.4) 


where 


-  r  cos  $ 


By  letting 


Equation  (C.1.4)  can  be  written  as 
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(C.1.5) 


(C.1.6) 


2c^jkR  +  ~ 


!fi_  i^i  +  i3c 

3s  ar2  r  dr 


(C. 1,7) 


Ue  now  use  separation  of  variables  to  solve  Equation  (C.1,7),  By 


substituting 


Cl(s.r)  -  S(s)X(r)  , 


(C.1.8) 


Equation  (C.1.7)  may  be  written  as 


2jkR  -  sr  -  ^  J  f  -  -l 


(C..1.9) 


where  is  an  unknown  constant. 


Therefore , 


v  >  2 

X'  '  +  tJ—  .  k  x  -  0 

r  1 


(C.1.10) 


and 
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f1  -  k*  -  2jkR 


(C. 1,11) 


Equation  (C.1.10)  may  be  re-written  as 


2  d  X  dX  2  2V  . 

:  —  *rS.«rX.O 

dr 


(C.1.12) 


and  the  solution  for  >  0  is 


X(r)  -  A1IQ(»c1r)  +  A2KQ(/c1r) 


(C.1.13) 


where  Aj_  and  A2  are  the  constants  ana  K0  and  Iq  are  the  modified  Bessel's 

function  of  zeroth  order  of  first  and  second  kind  respectively. 

2 

Note  1 ~  Negative  values  of  yield  functions  Jq  and  Yq  (Bessel's 

function  of  zeroth  order  of  first  and  second  kind  respectively)  which  do 

r/t  satisfy  the  boundary  conditions.- 
2 

Note  2 :  -  0  yields  a  logarithmic  function  which  again  does  not  satisfy 

the  boundary  conditions. 

For  large  distances  from  the  origin,  the  function  ci(s,r)  should  var,:.'»n 
Hence,  we  set  A]_  -  0  in  Equation  (C.1..13).  Therefore, 

X(r)  -  A,K  (k  r)  .  (C.1.14) 

v>  L 


The  solution  of  Equation  (C.1,11)  is  easily  written  as 


S(s)  -  Be 


(**-  2jkR)s 


(C. 1,15) 
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where  B  is  an  unknown  constant.  Hence  irons  Equation  (C.l.F),, 


C^s.r)  -  B^ 


(*^-2jkR)s 


Ko(<iD 


(C.1.16) 


where  B^  -  BA?  is  a  new  unknown  constant.  From  Equation  (C.!,2), 


2 

D  „  .  ,  KlS  i(k+<p-2kRs) 

c(s,r,t)  -  B..K  (K.v)e  eJ 
1  o  I 


(C.1.17) 


The  constants  B^  and  *3  evaluated  usirg  the  boundary  conditions  ere  found 
to  be 


B1  -  2,^  *nd  *1  -  "  ' 


Hence 


i(x.y.t)  -  K0(R  l/TTyW™-" 


(C.1.18) 
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APPENDIX  C . 2 

To  obtain  the  inverse  Laplace  Transform  of  the  function 


F(s.t) 


_ Qh  _±  jk _ 

s[j'2R  -jk  +  7s  +  jk 


(C.2.1) 


Let 


a  -  jk 


(C.2.2) 


b  -  72R  -  jk 


(C.2.3) 


F(s , t)  - 


0  Js  +  a 


s[7s  +  a  +  b] 


(C.2.4) 


F(s  ,  t)  -  Q(- 


7s  +  a  +  b 


s[Js  +  a  +  b]  s[7s  +  a  +  b) 


(C.2.5) 


F(s.t')  -  Q [ 


1 _ b_ 


s[7s  +  a  +  b] 


'I 


(C.2.6) 


Let 


H(s ) 


_ b _ 

s[7c  +  a  +  b] 


(C.2.7) 


or 


H(s) 


bf7s  -v  a  -  bl 
s(s  +  a  -  b^j 


20 


31  July  1986 
KR Ihz 


Denote 


a  —  a  -  b* 


(C.2.8) 


Then, 


H(s) 


b/s  ±  a 


s(s  +  a)  s(s  +  a) 


(C.2.9) 


uenote 


Z(s)  - 


b  Js  + 


S(5  +  Q) 


(C.2.10) 


Zisl 

b 


s  +  a  4-  b 


s(s  +  a)Js  +  a 


mi 

b 


s Js  +  a  s(s  +  a)Js  +  a 


ilsi 

b 


S Js  +  a  Q  s/s  +  a  (s  +  a)/s  +  a 


(C.2.11) 


Substituting  Equation  (C.2.11)  in  (C.2.9),  one  has 


Hvs)  - 


s  Js  +  a  asjs  + 


(s  +  a)J s  + 


s ( s  +  a) 


(C.2.12) 


Substituting  Equation  (C.2.12)  in  Equation  (C.2.6),  one  gets 
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F(s)  -  Q[; 


sjs  + 


asjs  +  a  a(s  +  a)J$  + 


s(s  +  a) 


(C.2.13) 


Taking  the  inverse  Laplace  Transform  on  both  sides,  one  finds 


f(x)  -  Q[~  -  ~  e  QXerfcbvx  -  erf^ax] 
a  a  a 


(C.2.14) 
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APPENDIX  C . 3 
Integral  evaluation 


Equation  (4.25)  is 

Ut) 

jjf  -  /  f(x,t)dt  ,  (C.3.1) 

0 


where 

f(x,t)  -  ^  (a  -  b^e  aXerfcb,/x  -  bja.  erf,/ax) 

Q 


(C.3.2) 


Substituting  Equation  (C.3.2)  in  (C.2,1), 

l(t)  £(t)  _  4(0 

~  ^  [  /  adx  -  b ^  /  e  QXerfcb,/x  dx  -  bja.  f  erfjax  dx] 

dt  a 

0  0  0 


(C.3.3) 


Let 


lit) 

1^  -  /  e  QXerfcb7x  dx 

0 


(C.3.4) 


and 


Ut) 

1 2  "  /  erf^ax  dx 
0 


(C.3.5) 


The  integral  in  Equation  (C.3.4)  may  be  evaluated  by  parts.. 
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r  rtvi?(t)  h2 

_  _  ,  e"aX  ^  1  -ax.  b  e ' D 

I,  -  erfcbyx  -  +  -  e  (-1)  ~  “ 

1  -  a  a  ,t 

L  J0  o  J*  J* 


T  1  e  erfcb7&  b  ,  e 

I  _  .  -  .  -  /  —  dx 

la  a  /  n  / 

a Jk  0  Jx 


(C.3.6) 


By  setting  ax  -  n2  in  Equation  (C.3.6),  we  can  show  that 


I  -  1  ■  erf M  , 

Lcr  a  r 

aja 


(C.3.7) 


The  integral  in  Equation  (C.3.5)  can  be  obtained  by  parts, 


-  x  evfjax 


i(t)  £(t)  _ 

/  “  cLA 

J  a  e 

*  I  x  _  — —  dx 
n  n  J*  Jx 


i.e.  , 


_  *(t) 


-  lerfJaZ  /  7> 


~  -  ax  . 

x  e  dx 


(C. 3.8) 


The  integral  in  Equation  (C.3.8)  may  be  evaluated  again  by  parts  and 
the  result  is 


lal  +V—  e  a^-  ~  erfjal 


I.  -  lerfjal  +V —  e  -  ~  evfjal 
2  if  a  2a 


(C.3..9) 


Hence  from  Equation  (C.3.3), 
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x  =  0  CAVITY  x  =  l(t) 


Figure  1.  Boundary  Conditions  for  Unsteady  Diffusion  Analysis. 
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